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ABSTRACT 

We present the results of completeness simulations for the detection of point sources as well 
as redshifted elliptical and spiral galaxies in the iiT'-band images of the Munich Near-Infrared 
Cluster Survey (MUNICS). The main focus of this work is to quantify the selection effects in- 
troduced by threshold-based object detection algorithms used in deep imaging surveys. There- 
fore, we simulate objects obeying the well-known scaling relations between effective radius 
and central surface brightness, both for de Vaucouleurs and exponential profiles. The results 
of these simulations, while presented for the MUNICS project, are applicable in a much wider 
context to deep optical and near-infrared selected samples. We investigate the detection prob- 
ability as well as the reliability for recovering the true total magnitude with Kron-like (adap- 
tive) aperture photometry. The results are compared to the predictions of the visibility theory 
of Disney and Phillipps in terms of the detection rate and the lost-light fraction. Addition- 
ally, the effects attributable to seeing are explored. The results show a bias against detecting 
high-redshifted massive elliptical galaxies in comparison to disk galaxies with exponential 
profiles, and that the measurements of the total magnitudes for intrinsically bright elliptical 
galaxies are systematically too faint. Disk galaxies, in contrast, show no significant offset in 
the magnitude measurement of luminous objects. Finally we present an analytic formula to 
predict the completeness of point-sources using only basic image parameters. 

Key words: surveys - infrared: galaxies - galaxies: photometry - cosmology: observations 
- techniques: image processing 



1 INTRODUCTION 

The luminosity function as well as the mass function, and to a lesser 
degree the number counts of galaxies provide an important obser- 
vational toolset for understanding the evolution of galaxies. Due to 
their statistical nature, these methods rely on the understanding of 
sample selection effects, i.e. the knowledge of what fraction of the 
true number of galaxies, as a function of their intrinsic properties, 
is actually present in the sample. In order to extract this information 
from the data, two different approaches are commonly used in the 
literature, creation of artificial objects or modification of observed 
objects. 



known, predictions about the total completeness of the data cannot 
be made. 



In the second ca se, a s conducted by e.g. Bershady et al. (1998 ), 



Arnouts et al. (1999) or Saracco et al. (199^), images of the ob- 



The synthetic objects are, as described by e.g. Martini (2001 ), 



usually point-sources or objects with galactic profiles, that are in- 
serted into the observed images. Then the fraction of objects recov- 
ered by the applied detection algorithm as a function of, e.g, the 
assigned apparent magnitude is computed. This approach exhibits 
two major drawbacks: The galaxies are created using a discrete set 
of half-light radii and a continuous range of magnitudes, thus ignor- 
ing known relations between su rface-brightness aii d effective ra- 
dius, Jikedie^undamental_^^ et al. 
2000). The results show the completeness limits for distinct types 
of objects, but as the actually observed ratio of these types is un- 



served objects are dimmed or brightened to produce artificial ob- 
jects. The advantage of this approach is, that it preserves the ob- 
served mix of galaxies, assuming that it remains unchanged with 
time, yet the resulting completeness fractions for extended objects 
are questionable, as long as the resolution of images is seeing- 
limited. When the bright extended objects are dimmed, they are still 
profile dominated, whereas in reality a faint object would be see- 
ing dominated. The same holds true for artificially brightened faint 
sources. The resulting objects would still be seeing dominated. Fur- 
thermore, the observed size of a local object of given magnitude is 
different from the size of a distant object of the same brightness. 
Consequently this method only yields information about the prob- 
ability to detect an object if it was fainter, but not about the proba- 
bility to detect an intrinsically faint object. 

Another drawback of both approaches is that the obtained re- 
sults can only be used to define the completeness of the survey in 
apparent magnitudes. The magnitudes and radii of the artificially 
constructed objects all are physically plausible, but the objects do 
not occupy one plane in the (/ie) — re — jz parameter space. There- 
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fore the resuUs derived using galactic profiles are not usable to cor- 
rect absolute magnitudes, as effects of the completeness of different 
galaxy types at given redshifts are unknown. 

In deep extra-galactic surveys, the observed galaxies span a 
wide range in intrinsic profile shape, intrinsic brightness, and intrin- 
sic size. The apparent quantities vary with cosmological distance, 
such that the fraction of galaxies visible is also a function of red- 
shift. The need to simulate objects obeying the known scaling rela- 
tions for galaxies was pointed out several times in the literature. Us - 
ing profiles obeying the magnitude-radius relations, ibshii (1993) 
predicted the number of objects lost in number count analysis. He 
finds a strong dependence of the detection rate on the applied de- 
tection criteria and magnitude measurement algorithm, leading to a 
larger number of undetected faint galaxies at high redshifts. 

Dalcanton (1998) analysed the biases of the luminosity func- 
tion introduced by the cosmological effects of size variation with 
redshift, and cosmological dimming for galactic profiles in depen- 
dence of size and magnitude, taking into account effects of seeing. 
Starting from the deficiency that magnitudes are usually measured 
as some sort of isophotal magnitude, that is directly influenced by 
the above mentioned effects, she finds the possibility of a severe 
underestimation of the true luminosity function introduced by the 
fact that an object and distance dependent part of the light is lost 
outside the limiting isophote. 

The simulations presented in this paper were performed by 
adding artificially created objects into the i^'-band images of 
the Municl^^Jear-Infrare^^^^ et al. 

2001). We analysed the detection probability and the reliability 
for recovering the assigned magnitude for point-like sources (Mof- 
fat profiles), elliptical galaxies (de Vaucouleurs profiles), and spi- 
ral galaxies (exponential profiles). The radii and magnitudes of 
the simulated galaxies were distributed accordi ng to the projected 
fundamental-plane relation ( Bender et al. 1992 ) for ellipticals and 
a Freeman law type relation ( Freemairi97(]| ) constructed from ob- 
served data for spirals. The galaxies were simulated at five dis- 
tinct redshifts between z = 0.5 and z = 1.5, taking into account 
size variation with redshift as well as cosmological dimming and 
K-correction. A flat universe with Q.m — 0.3, JIa = 0.7 and 
Ho = 65kms"^ Mpc"^ was assumed. 

The results of the completeness simulations for the different 
fields of the MUNICS survey presented here, illustrate the princi- 
pal limitations of imaging surveys. A full correction of the incom- 
pleteness would only be possible if the type mix of galaxies was 
known. The presentation of the results here is limited to the K'- 
band data, but as the selection biases are caused by the physical 
nature of the objects, they are applicable in a much wider context 
to deep extragalactic surveys spanning the optical and near-infrared 
wavelength regimes. The conclusions drawn here result from the 
use of a threshold-based detection algorithm. Accordingly they will 
hold true for other similarly created datasets as well. 

Finally, a reliable and handy analytic formula to estimate the 
completeness limit of a survey for point-like sources is presented. 
In Section ^ the implementation of the simulations and the gen- 
eration of the artificial objects is described. Section ^ presents the 
results of the simulations, and a discussion of these. To compare the 
results with analytic predictions, we pre sent the results of an analy - 
sis using the visibility theory devised by Disney & Phillipps (1983) 
in Section ^ for the detection probability and in Section ^ for the 
lost-light fraction. In Section pi we discuss the effects of seeing in 
the above mentioned analysis. Finally, in Section ^ we present the 
results of the completeness simulations for the MUNICS fields. 



1.1 The Munich Near-Infrared Cluster Survey (MUNICS) 

The Munich Near-Infrared Cluster Survey (MUNICS) is a wide- 
field medium-deep survey in the near-infrared K' and J pass-bands 
(Drory et al. 2001; MUNICS 1 hereafter). The survey consists of a 
./^'-selected catalogue down to K' < 19.5 covering an area of 
1 square degree. Additionally, 0.35 square degrees have been ob- 
served in /, R, and V. The layout of the survey, the observations 
and data reduction are described in Drory et al. (2001 ). 



2 IMPLEMENTATION OF THE SIMULATIONS 

Detection probabilities and photometry results were analysed for 
three different profile shapes: de Vaucouleurs profiles, exponential 
disks and point-like sources. 

For each profile type and each image, 200 artificially created 
objects were added to the image, taking into account the noise 
properties of the background and the photometric zero-point. The 
resulting image was processed in the same way as described in 
MUNICS 1 regar ding detectio n and photometry, using the software 
package YODA (Drory 2002). The resulting object catalogue was 
used to calculate the fraction of recovered objects and to com- 
pare the resulting photometry in Rron-like elliptical apertures to 
the quantities in the input catalogue. 

This procedure was repeated 500 times resulting in a total of 
100 000 artificial objects per type. In each run, the artificial ob- 
jects were distributed randomly in (x,y)-position across the image, 
excluding a 25 pixel wide strip along the image borders, and re- 
quiring 20 pixel distance to existing objects, as we intentionally 
avoid crowding which would introduce effects beyond the scope of 
this work. Extended sources were convolved with a Gaussian with 
a FWHM of the measured seeing in the image. Poisson noise was 
added to all profiles. In the following we discuss the generation of 
objects of each of the examined profile types in more detail. 



2.1 Stars and point-sources 

Poin t-like sources were simulated using a Moffat profile (M offat 
1969) of the form 

Ic 



7(r) = 



(1) 



with the characteristic intensity Ic- The size parameter R, defin- 
ing the radius of the created object, was set to the measured value 
from point-like real sources in the analysed image, and the Moffat 
parameter /3 was fixed at the canonical value of 2.5. The apparent 
magnitudes assigned to the point-sources were chosen randomly 
from a constant probability density in the range 15 < mj^i < 25. 



2.2 Elliptical galaxies 

Elliptical galaxies wer e simulated using a de Vaucouleurs profile 
(de Vaucouleurs 1948). For a galaxy with effective radius re and 



effective intensity Je (defined as the intensity at re), this profile can 
be written as 



/(r) = Jeexp -7.67 



(2) 



For the sake of simplicity an r^^* law is applied, ignoring 
relations between luminosity and the shape of the profile. But it 
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de Vaucouleurs 


exponential 


0.5 


-0.63 


-0.68 


0.75 


-0.73 


-0.83 


1.0 


-0.84 


-0.97 


1.25 


-0.88 


-1.03 


1.5 


-0.87 


-1.02 



Table 1. K-coiTections in magnitudes for de Vaucouleurs and exponential 
profile s applied in th e simulations. The values were derived from empirical 
SEDs (iBender 200ll). 



should be kept in mind, that an r^^" profile with n > 4 would 
show the same biases as described here. To create a realistic popu- 
lation of elliptical galaxies, the radii and absolute magnitudes of the 
galaxies were distributed according to the Kormendy relation (Ko- 
rmendy 1985), a proj ection of the local fundamental-plane relation 
( Burstein et al. 199 7|i. We us ed the A'-band Kormendy relation pub- 
lished by Pahre et al. (1995 ) for the effective radius of the galaxy 
i?e in kpc and the mean effective surface brightness (/le) in mag 



logTJe = 0.332{^ie) - 5.090. 

Using the definition of (/le) 

m = —1.995 — 5 log re + (/ie). 



(3) 



(4) 



equation (|3|) can be transformed into a relation between absolute 
magnitude M and effective radius Re. 



-1.99 log i?e - 23.235. 



(5) 



The absolute magnitudes were chosen randomly and uniformly in 
the range —26.5 < Mx' < —20.5, corresponding to a range of ±3 
magnitudes aro und for the local K-band luminosity function 
dLoveday 200C|). 



To simulate a galaxy at redshift z, the apparent radius was 
calculated from the given physical radius using the angular dis- 
tance, and the absolute magnitude was transformed into the ap- 
parent magnitude using the luminosity distance. During each rep- 
etition of the simulation 40 galaxies for each of the five redshifts 
z £ {0.5, 0.75, 1.0, 1.25, 1.5} were created. The axis ratios e were 
selected randomly in the range 0.7 < e < 1.0, the position angles 
were chosen arbitrarily. For the galaxies ai z = 0.5 the absolute 
magnitude range was shifted —25.5 < M^' < —19.5 to better 
trace the faint-end drop-off of the completeness curve. The result- 
ing apparent magnitudes were a djusted using o nly K-corrections 
derived from model SEDs from Bender (2001). The used values 
are listed in Table [l[ Possible effects introduced by luminosity evo- 
lution with redshift were explored, and are discussed at the end of 
section H. 



2.3 Spiral galaxies 

Spiral galaxies were simulated as pure exponential profiles with 
no bulge component. An anal ysis of the bulge-disk decomposition 
for the galaxies published by Gavazzi et al. (200C ) shows that the 
contribution of the bulge to the total luminosity is less than 18%. 
We therefore neglect the contribution of the bulge. 

For a galaxy with half-light radius Vh and central intensity le 
the intensity profile can be written as 



I(r) = Joexp f-1.6783 — 



(6) 
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Figure 1. Absolute ii''-band Magnitudes A/^ and effective radii for the 
galax ies fitted bv a pure exponential profile published by G avazzi et al. 
(2000) (filled squares) and some during the simulations created artificial 
objects (open triangles). 



galaxies fitted by a pure exponential profile from ff-band sur- 
face p_hotoinetr^_data^fs£ira^^ala^^ et al. 
(2000) were used. The absolute J?-band magnitudes were trans- 
formed into the if' -band using H — K ' colours derived from em- 
pirical SEDs presented in Bender (2001 ). The mean correction was 
H — K' ~ 0.128 mag. Using the objects' distances, the half-light 
radii in kpc Rh were computed from the published values. The re- 
sulting distribution of galaxies in the log Rh — Mj^i plane was then 
approximated by assuming 



-5 (log i?h ±0.25) - 19.7. 



(7) 



For a population of local spiral galaxies, this distribution would be 
analogous to a Freeman law with a AT' -band central surface bright- 
ness of 17.52 mag arcsec"^ and an rms of 1.73 mag. 

Fig. |l| shows the absolute magnitudes and effective radii of 
the Gavazzi et al. data and of some artificially created objects. The 
distribution of the artificial objects reproduces that of the observed 
population reasonably well, extending to somewhat brighter mag- 
nitudes. 

The absolute magnitudes of the objects were randomly cho- 
sen in the range —26.5 < M^' < —20.5, again corresponding 
to a range of ±3 magnitudes ar ound Af^ for a local K-band lumi- 
nosity function (Loveday 200C). Using the same techniques as for 



To create objects with realistic magnitude-size ratios, all 



the elliptical galaxies, the spirals were simulated for the same five 
redshifts, with axis ratios in the range < e < 1 and arbitrary 
position angles. The resulting apparent magnitudes were adjusted 
using onl y K-corrections derived from model SEDs shown in B en- 
der (2001), no evolution corrections were applied. The results of 
tests for influences of luminosity evolution are discussed at the end 
of section ^ The applied values are listed in Table |l| The mean 
surface brightnesses (/i) were derived from a sample of local spiral 
galaxies. Brightness evolution with redshift would make the galax- 
ies brighter, while the radii are known not to change significantly 
up to redshifts of z ~ 1. 
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rriK' [mag] 

Figure 2. Probability Pj:){m) to re-detect an artificially created point- 
source (upper panel), an extended object with a de Vaucouleurs profile (mid- 
dle panel) or exponential profile (lower panel) as a function of the objects' 
K'-band input magnitude mjfi . The different line types in the middle and 
lower panel show the completeness for the five analysed redshifts: z = 0.5 
(soUd), 0.75 (dash), 1.0 (dot), 1.25 (dashdot) and 1.5 (dashdotdot). 



rriK' [mag] 

Figure 3. Mean magnitude difference 5m = mj^i — rrimeas between the 
assigned magnitude / of an object and the Kron-magnitude measured 
by the photometry software rrimeas in bins of 0.25 mag as a function of the 
input i^'-band magnitude mjf/ . The top panel shows the results for point- 
sources, the middle panel for de Vaucouleurs profiles and the lower panel 
for exponential profiles. The errorbars show the standard deviation of 5m, 



3 RESULTS AND DISCUSSION 

The results for one of the MUNICS fields with a seeing of approx- 
imately one arcsec FWHM are shown in Fig. ^. The upper panel 
shows the results for the point-sources, the middle and the lower 
panel for the de Vaucouleurs and exponential profiles, respectively, 
for the five redshifts simulated. An interesting effect seen in the fig- 
ure is that, for higher redshifts z>l, the detection probability does 
not reach unity even for the brightest objects, forming a plateau at 
some lower value. This effect is mainly caused by the cosmologi- 
cal surface brightness dimming and discussed in detail in Section ^ 
below. Fig. |^ shows that the completeness for point-like sources 
provides a rough but reasonable approximation for the complete- 
ness of the analysed extended objects up to a redshift of 2; = 1.0. 

For each detected object, the difference between the in- 
put magnitude and the measured Kron-magnitude was computed. 
Fig. ^ shows the mean magnitude differences for the analysed pro- 
files, averaged in bins of 0.25 mag and with the standard deviation 
of the measurements in the bin indicated as errorbars. 

The object recovery probabilities for the high-redshift de Vau- 
couleurs profiles, and to a lesser degree for the exponential profiles 
exhibit a significant detection bias compared to lower redshift ob- 
jects, as shown in Fig. ^. The fact that the objects at these redshifts 
never reach a detection probability of one is caused by the distribu- 
tion of their physical parameters along the fundamental-plane rela- 
tion - in case of the ellipticals - and according to the Freeman law 
- in case of the disks. In both cases, as the object's size increases 



its radius grows and therefore the average surface brightness de- 
creases. As a result, even the brightest objects of the distribution 
fail to produce a large enough area above the threshold in surface 
brightness to be detected with high probability in the presence of 
noise. In Section ^ these results are compared with the theoretical 
predictions of the visibility theory, which will confirm the above 
statement. 

The detection deficiencies at low redshifts found for the ex- 
ponential profiles are caused by the flat light distribution of the 
objects, resulting in a low central surface brightness even for the 
luminous objects. The image noise can then scatter these objects 
below the detection thr eshold. This result is i n general agreement 
with the predictions of McGaugh et al. (1995 ). In their pape r they 
re-derive the formalism devised by pisney & Phillipps (1983 ) using 
two parameters for the disk galaxies. They predict that disk galaxies 
would suffer from detection biases introduced by their low central 
surface brightness. 

iviartmi (ZUUl) performed simulations similar to the ones pre- 
sented here, analysing the detection probabilities of Moffat, de Vau- 
couleurs and exponential profiles. The galactic profiles were cre- 
ated for a discrete set of half-light radii and a continuous range of 
apparent magnitudes. The simulations show a decrease of the lim- 
iting magnitude with increasing half-light radius. This effect can 
be found in the simulations presented here, by comparing the com- 
pleteness magnitude at different redshifts. With increasing redshift, 
the sampled range in absolute magnitudes is shifted toward fainter 
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de Vaucouleurs profile exponential profile 




r/re 

Figure 4. Intensity integrated over image pixels for a de Vaucouleurs profile 
(solid line) and an exponential profile (dashed line) at redshift z = 0.5 with 
similar total magnitude and effective radius. The upper left panel shows the 
case in the absence of seeing, the upper right one, with a seeing of 0.5 arcsec 
and the lower left with a seeing of 1.6 arcsec. Plotted is the logarithm of the 
intensity against the radius relative to the effective radius. 



rriK' [mag] rriK' [mag] 

Figure 5. Results of a simulation assuming luminosity evolution with red- 
shift. Shown are the simulation results for de Vaucouleurs profiles (left col- 
umn) and exponential profiles (right column) at a redshift 2 = 1 for a seeing 
of 0.8" (upper row) and 1.6" (lower row). The solid curve shows the results 
obtained from the normal simulations, the dashed line for a simulation with 
objects brightened artificially by 0.75 mag. 



apparent magnitudes, resulting in an increase of effective radius 
at given apparent magnitude. The detection bias against liigh red- 
shifted elliptical galaxies we showed cannot be found by the simu- 
lations performed by Martini, as these effects are only visible when 
distributing the objects along a fundamental-plane relati on. 

The bias first shown by Disney & Phillipps (1983 ), that with 
increased distance more low-luminosity system would be lost from 
the survey, is reproduced by our simulations as well. With in- 
creasing redshift, the magnitude of the completeness limit stays 
approximately constant, or gets brighter, while the range of sam- 
pled absolute magnitudes is shifted to fainter apparent magnitudes. 
This means, that with increasing redshift, the observable range in 
absolute magnitudes moves toward higher luminosities, the low- 
luminosity systems become unobservable. 

The deviations of the measured magnitude from the true mag- 
nitude for the de Vaucouleurs profiles as seen in Fig. ^, can be ex- 
plained as resulting from the estimate of the Kron-like aperture ra- 
dius under the conditions of a surface brightness limited detection 
procedure and the involved intrinsic brightness profile of the ob- 
jects. The upper left panel of Fig. Q shows the intensity of a de 
Vaucouleurs and an exponential profile with similar total magni- 
tude and effective radius, as a function of the radius in units of the 
effective radius in absence of seeing. Assuming that both profiles 
are detected at a similar intensity level, the Kron-radius of the de 
Vaucouleurs profile would be smaller compared to the exponential 
profile, even though the effective radius of both objects is the same. 
The measured light within the underestimated radius results then in 
a too faint object magnitude. Combined with the previously men- 
tioned surface brightness distribution along the fundamental plane, 
this leads to an increased amount of lost flux and a larger error 
in the output magnitude for brighter objects, as these have lower 
surface brightnesses. In Section ra we will confirm this explanation 
with predictions based on computations of the visibility function 
and the lost-light fraction. 



Our results confirm the predictions of McGaugh (1994), that 
the magnitudes of distant spiral galaxies would be measured cor- 
rectly for luminous galaxies, and slightly underestimate them for 
the systems with low-luminosity. 

The deviations of the measured from the input magnitu de we 
find, are compatible with the simulations of Martini (2001), who 
compared the reliability of different photometric methods for point- 
sources and exponential profiles. For Kron-like magnitude mea- 
surements both shapes would be measured correctly at the bright 
end, and slightly underestimated at the faint end. 

To explore the effects caused by the increase of the objects' 
surface brightness introduced by luminosity evolution with red- 
shift, an additional set of simulations for objects at z = 1 was cre- 
ated. The brightnesses of the created artificial de Vaucouleurs and 
exponentia^^rofites^asjncreased^y0^7£^^ Propris 
et al. 1999). The results shown in Fig. |5| show no significant change 
of the magnitude of the 50% completeness limit. The reason for this 
is that the increase of the detectable area is small in the magnitude 
range of the 50% limit. The increase in size is larger for brighter 
magnitudes, but the positive effect is mainly lost due to the already 
high detection probability, or may result in a slightly higher plateau 
level. 



4 VISIBILITY THEORY 



Disney Sl Phillipps (1983) and Phillipps et al. (1990) analysed the 



dependence of the visibility of a galaxy on its surface brightness 
and apparent radius in a survey with given detection constraints in 
surface brightness and object radius. They estimated the maximum 
distance at which a galaxy with a given magnitude and effective 
radius can be seen, by calculating the distance at which the surface 
brightness at the limiting radius drops below the detection thresh- 
old. 

Even though our goal is not to predict the maximum distance 
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Figure 7. Predictions of the visibility function in comparison witli the results of completeness simulations for de Vaucouleurs profiles (upper row) and 
exponential profiles (lower row) at redshift 2 = 0.5 (left column) and z = 1.25 (right column). The upper panel of each plot depicts the area above detection 
threshold. The thin solid line shows the area above the detection threshold integrated over the image pixels in absence of seeing, the thick solid line the same 
but in presence of appropriate seeing. The smooth dotted curve gives the area above the surface brightness limit in absence of seeing, the dashed line in 
presence of seeing. The thin horizontal line is the area of the minimum number of connected pixels needed to detect an object. The lower panels show the 
object recovery probability as a function of assigned input magnitude measured for a MUNICS image with seeing of ~ 1 arcsec. 



out to which a given galaxy type can be seen, but to analyse biases 
that occur when observing galaxies at high redshifts, this theoret- 
ical approach can be used here as well to predict the behaviour of 
the used threshold-based detection algorithms. 

To detect an object we require a minimum number of consec- 
utive pixels to lie above a given brightness threshold, usually ex- 
pressed in units of the standard deviation of the background noise. 
The area required to be above the threshold is usually determined 
by the seeing disk (resolution element) size. Both values are ad- 
justed such that faint real objects are detected at a tolerable rate of 
false detections. 

In our case the minimum number A'^ of consecutive pixels is 
chosen to be A'' = 7r(1.4s/2)^, with s being the seeing FWHM in 
the image. The threshold t is set to three times the standard devi- 
ation a of the local background noise, for details see D rory et al. 
(2001). 

In the case of circularly symmetric profiles, the calculation of 



the area above t is trivial, and the limiting surface brightness ^nm 
may be written as 

= m^p ~2.^\og-^ (8) 

for the magnitude zero-point niz-p and the pixel scale p (0.396 arc- 
sec/pixel in MUNICS). 

To create comparable completeness curves from the simula- 
tions as discussed above, a set of simulations for point-sources with 
Moffat profiles and circular face-on galaxies were calculated. These 
results are shown in Fig. ^ for point-sources and in Fig. ^ - in the 
lower panel of each plot - in the upper row for de Vaucouleurs 
profiles, and in the lower row for exponential profiles at redshifts 
z = 0.5 andz = 1.25. 

For Moffat profiles, the radius R of the circular area A above 
the limiting isophote can be calculated as 



R^8\l 10^^T3r^ _ 1 (9) 
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Figure 6. Area above detection threshold for point-sources with Moffat pro- 
files. Theu^ESLE^IlsLMlSSliil'^ predictions derived from the visibility the- 
ory (Disney & Phillipps 1983). The thick solid histogram shows the area 



above the detection threshold integrated over the image pixels in presence 
of seeing. The smooth dashed line gives the area above the surface bright- 
ness limit as predicted by the inversion of the profile. The thin horizontal 
line is the area of the minimum number of connected pixels needed to de- 
tect an object. The lower panel shows the object recovery probability as a 
function of assigned input magnitude. 



with seeing s and characteristic surface brightness Hc- In the up- 
per panel of Fig. ^ the dashed hne shows the area A resulting from 
the analytic calculation following equation the thick solid his- 
togram gives the area integrated over the image pixels for a simu- 
lated object. The thin horizontal line indicates the minimum limit- 
ing area needed to detect an object. 

Equation (B) can be transformed into 



2.5 loB 



V 



(i+(f)T 

/3-1 



(10) 



Using equation 



this can be written as 

2p - 2.5 log 



tns- 



(i+(f)T 



v 



p2 (f3 - 1) 



(11) 



These formulas provide a simple way to estimate the completeness 
limit for point-like sources for a given image, using easily mea- 
surable parameters. Tests using the MUNICS data have shown that 
this formula provides a robust estimate of the ~ 50% completeness 
level. 

For de Vaucouleurs profiles with effective radius and effec- 
tive surface brightness fie within r^, 7? can be written as 



and as 



. 8.3268 



+ 1 



(12) 



R = rh 



MO 



822 



(13) 



for exponential profiles with half-light radius and the central 
surface brightness fio- 

The dotted lines in the upper panels of Fig. ^ show the area 
within the limiting isophote as a function of apparent magnitude, 
as calculated using equation ( p^ and Equation ( p^ for de Vau- 
couleurs profiles and exponential profiles, respectively. For com- 
parison, the same values extracted from objects in the simulations 
are plotted as solid lines. The area corresponding to the minimum 
area required for detection is indicated as a horizontal line. The in- 
tersection of this line with the curve then provides an estimate of 
the completeness limit. 

The lower panels of these figures show the detection proba- 
bility as a function of apparent magnitude. Results are shown both 
with and without seeing. Seeing was modelled as a convolution 
with a Gaussian profile. In the case of the de Vaucouleurs pro- 
files, moderate seeing ( ^ 1 arcsec) improves detectability, since it 
distributes flux outwards, converting the steep core of the de Vau- 
couleurs profile to a larger and flatter flux distribution. This effect 
is much weaker for disk profiles since these are flatter, anyway. 
Note that worse seeing again worsens the detectability of objects, 
since at some value, it will distribute too much flux outwards and 
cause the core of the profile to drop beyond the threshold. For a full 
discussion, see Section ^. 

From these plots it becomes apparent, that the theoretical pre- 
dictions provide a good estimate of the completeness function (as 
the detectable area drops below the required value at the same ap- 
parent magnitude at which the completeness function drops off), 
but only provided that seeing is taken into account. 

In the case of the high-redshift {z — 1.25) galaxies, the detec- 
tion probability never reaches one. This results from the fact that 
even for the brightest simulated objects ~ with an absolute mag- 
nitude 3 mag above A'l* , but lower mean surface brightness ~ the 
object's core lies barely above the detection threshold for all mag- 
nitudes. Therefore, noise easily scatters objects below the threshold 
and the detection probability is always smaller than unity. 

Looking at the predictions of the visibility function for the ex- 
ponential profiles at high redshifts, the low detection probability 
surprises. Although the detectable area is much larger than the lim- 
iting one, the completeness fraction remains low. Due to its rather 
flat light-distribution, the exponential profile is much more suscep- 
tible to distortions caused by noise than the steeper de Vaucouleurs 
profile. These scattered noise pixels can then either reduce the ob- 
jects' size below the limiting radius, or form additional false max- 
ima leading to a mis-detection in the form of several objects, not 
recognising the artificially created one. 



5 LOST-LIGHT FRACTION 

The magnitude differences between the input and the measured 
magnitudes shown in Fig. ^ exhibit a strong deviation for ellipti- 
cal galaxies at the bright end of the distribution. The results from 
the analysis of the visibility function can be used to calculate the 
lost-light fraction. Assuming that an object is detected out to the 
limiting radius rnm where the surface brightness drop s below the 
detection limit, the intensity weighted Kron radius Vk ( Kjon 1980) 
can be calculated for circularly symmetric objects as 



(14) 
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Figure 8. Deviation of tlie measured magnitude from tlie assigned input magnitude 5m for de Vaucouleurs profiles (left figure) and exponential profiles (right 
figure) with a seeing of 0.8 arcsec, as a function of the assigned _fC'-band input magnitude. To provide a clue for the width of the distribution, the rms is shown 
as errorbars. The solid line shows the theoretical prediction of the lost-light fraction in absence of seeing, the dashed line for profiles convolved with a Gaussian 
seeing. The horizontal scale (dashed-dotted line) at 5m = 1 in each panel gives the effective radii of the simulated objects for Re £ {10, 5, 2, 1, 0.5, 0.1} kpc 
for the de Vaucouleurs profiles at 2 = 0.5, re G {40, 20, 10, 5,1,0.1} kpc for the other elliptical galaxies and re G {20, 10, 5,2,1} kpc for the exponential 
profiles. The large filled triangle denotes the apparent magnitude and radius of an M* galaxy. 



The galactic profile is then integrated out to some factor (in 
our case 2.5) times the Kron radius, and the corresponding mag- 
nitude is calculated. Fig. ^ shows the results of this calculation of 
the lost-light fraction for circular symmetric de Vaucouleurs and 
exponential profiles in comparison with the results obtained from 
the completeness simulations. 

The measured magnitudes of the elliptical galaxies deviate 
strongly from the assigned input magnitudes due to underestima- 
tion of the Kron radius caused by the steep decline of the profile 
in the central detectable part. In case of the exponential profiles 
the differences in the photometry are much smaller, since the flat- 
ter profile causes the Kron-radius estimate to be closer to the true 
value, even if only the inner part of the profile is detected. 

The analytically calculated lost-light fraction predicts a 
slightly lower difference than the measurements show (except for 
the highest-redshift bin). This can be explained by the fact, that the 
theoretical approach is based on the ideal case, where the object 
is detected exactly out to the maximum visible radius, and that the 
form of the profile is not distorted. In reality, the area is pixelated 
and integration over the pixels in calculating the Kron radius causes 
smaller values. Additionally, noise causes mis-measurement of the 
object size and shape. 

In contrast, in the highest redshift bin (z — 1.5) analysed here 
in the context of MUNICS, the analytically calculated lost-light 
fraction predicts higher deviations than actually measured in the 



simulated images. These objects are only detectable due to addi- 
tional pixels b eing scattered above the threshold by noise (as can 
be seen in Fig. hol which shows that the whole profile is below the 
threshold). This additional flux causes the photometry to yield too 
bright magnitudes. 

The same effect causes the photometry to yield too bright 
magnitudes for objects at the faint end of the luminosity function 
in all redshift bins. These are also too faint to be detected without 
the presence of pixels scattered above the threshold by noise. 

As the mean surface brightness increases towards fainter ob- 
jects, and as those are also intrinsically smaller, a larger fraction of 
the total profile is visible around M* , and the measurement of the 
total magnitude improves. 

It should be kept in mind that the largest differences occur 
for the rather rare objects with absolute magnitudes 3 mag brighter 
then AI* , while the magnitudes of M* objects are measured cor- 
rectly. In the calculation of the luminosity function, this would 
cause the brightest objects to be redistributed to lower absolute 
magnitudes. But again, the numbers of such bright galaxies are 
rather low (~ 10~* times lower space density than M* objects), 
and thus the contamination of the LF by these galaxies should be 
rather small. At the faint end, the situation is similar Although 
these objects are much more numerous, they are actually below the 
detection limit. Therefore these are only detected by chance (about 
200 out of 20000 simulated objects at z = 0.5 were detected be- 
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Figure 9. Deviation of tlie measured magnitude from tlie assigned input magnitude Sm for de Vaucouleurs profiles (left figure) and exponential profiles (right 
figure) with a seeing of 1.6 arcsec. To provide a clue for the width of the distribution, the rms is shown as errorbars. The solid line shows the theoretical 
prediction of the lost-light fraction in absence of seeing, the dashed line for profiles convolved with a Gaussian seeing. The horizontal scale (dash-dotted line) 
at 5m = 1 in each panel gives the effective radii of the simulated objects for Re £ {10, 5, 2, 1, 0.5, 0.1} kpc for the de Vaucouleurs profiles at z = 0.5, 
re G {40, 20, 10, 5, 1, 0.1} kpc for the other elliptical galaxies and 6 {20, 10, 5, 2, 1} kpc for the exponential profiles. The large filled triangle denotes 
the apparent magnitude and radius of an M* galaxy. 



yond 20 mag), Thus the influence of their wrongly measured mag- 
nitude on statistics like the luminosity function is almost negligible. 



6 THE EFFECT OF SEEING 

To explore the influence of seeing on the visibility, seeing con- 
volved galactic profiles were used to calculate the visibility func- 
tion and the lost-light fraction. The seeing was simulated using a 
two-dimensional Gaussian of the form 



/(r) 



27r(T" 



■ exp 



r 

2^ 



(15) 



The width of the Gaussian kernel a was calculated from the width 
s of the measured seeing PSF as 

(16) 



\/8Tn2 



The galactic profiles Ig were convolved with the Gaussian resulting 
in the convolved profile I'g 



oo 27r 

I» = J J r'h{r)e 





dipdr' 



(17) 



Using convolved profiles, the calculation of the visibility func- 
tion and the lost-light fraction were repeated on images with two 
different seeing values of 0.8 and 1.6 arcsec. 

Fig. [l^ shows the results of the completeness simulations and 
the calculation of the visibility function for de Vaucouleurs profiles 
with seeing of 0.8 and 1.6 arcsec, Fig. [ill the same for exponential 
profiles. 

An increase of the seeing distributes more light outward from 
the central parts of the profile within the seeing disk, resulting in 
a smoothed light distribution in the centre, as shown in Fig. ^. As 
discussed above, moderate seeing improves detectability of faint 
objects by distributing light from the bright centre more evenly 
across a larger area without reducing the light in the centre below 
the threshold. As the seeing gets larger, this effect is counter bal- 
anced by the fact that now the amount of light redistributed away 
from the centre becomes so large that the central parts of the profile 
falls below the detection threshold. 

Fig. |l^ shows this effect for the case of the steep de Vau- 
couleurs profile The detection limit is fainter by roughly one mag- 
nitude in the presence of a seeing of 0.8 arcsec compared to the 
seeing-free case. On the other hand, once the seeing is as large as 
1.6 arcsec, the detection limit has dropped close to the no seeing 
case. In the highest-redshift bins, detectability of elliptical galax- 
ies is depressed below the no-seeing case, since now the seeing- 
convolved profile is entirely below the detection threshold. 
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Figure 10. Visibility function and completeness function for circularly symmetric de Vaucouleurs profiles convolved with a Gaussian seeing of 0.8 arcsec (left 
figure) and 1.6 arcsec (right figure) at redshifts z £ {0.5, 0.75, 1.0, 1.25, 1.5}. The upper panels show the theoretical predictions of the visibility function in 
the analytic case (dashed) and integrated over the image pixels (solid). The lower panels show the object recovery fraction as a function of input magnitude. 



The exponential profiles - shown in Fig. [H] - suffer from the 
same effects, resulting in decrease of the completeness limit with 
increased seeing. The low overall recovery fraction of exponen- 
tial profiles compared to de Vaucouleurs profiles at high redshift is 
caused by the flatter light distribution and the lack of a central peak 
compared with the de Vaucouleurs profiles. Due to this, the ob- 
jects' appearance is much more irregular in the images due to the 
significant fraction of flux in poissonian noise. This results, on the 
one hand, in pixels dropping below the detection limit, and, on the 
other hand, false maxima, resulting in the detection of two or more 
(distinct) structures. These effects have the strongest impact for the 
circularly symmetric (face-on) objects used here. In the more real- 
istic approach using randomly distributed ellipticities - as shown in 
Fig. g - the impact of these effects is much less significant, due to 
the steeper apparent profiles for inclined disks. 

The lost-light fraction - shown in Figs. |^ and |^ - for objects 
with seeing, in contrast, shows no significant change with the in- 



crease of the seeing except close to the completeness limit and at 
high redshift. 



7 K' -BAND COMPLETENESS LIMITS FOR THE 
MUNICS SURVEY 

Here we present the results of the c ompleteness simula tions for the 
fields of the MUNICS survey (see Drory et al. 2001 ) in the K'- 
band. To parameterise the shape of the completeness curves a com- 
bination of two power-laws with the functional form 



Poim) 



Po 



(18) 



+ 



was fitted to the results. The best-fitting parameters pp, rno, a and 
/3 are given in Table ^ for point-like sources. Tables^ and ^ for de 
Vaucouleurs profiles and Tables q and H for exponential profiles. In 
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Figure 11. Visibility function and completeness function for circularly symmetric exponential profiles convolved with a Gaussian seeing of 0.8 arcsec (left 
figure) and 1.6 arcsec (right figure) at redshifts z £ {0.5, 0.75, 1.0, 1.25, 1.5}. The upper panels show the theoretical predictions of the visibility function in 
the analytic case (dashed) and integrated over the image pixels (solid). The lower panels show the object recovery fraction as a function of input magnitude. 



case of the latter ones the parameters are given separately for the 
five analysed redshifts. 



8 SUMMARY 

We presented the results of extensive completeness simulations for 
imaging surveys of faint galaxies, taking into account their known 
size - surface-brightness relations. The absolute magnitudes of the 
elliptical galaxies simulated as de Vaucouleurs profiles following a 
local fundamental-plane relation. The disk galaxies represented by 
pure exponential profiles were modelled according to a Freeman 
law. These physical parameters were converted into apparent sizes, 
depending on the simulated redshift of the object. The simulations 
were carried out for the K'-hwA images of the MUNICS survey, 
but the deficiencies found should remain valid in other comparable 
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Table 2. Parameters of the best fit of equation (|l8|) to the results obtained 
from the completeness simulations for point-like sources. The normalisa- 
tion PO provides a rough estimate of the height of the plateau the detection 
probability reaches. The break magnitude mo can be used as an approxi- 
mation of the 50% completeness limit, if a higher maximum detection rate 
is reached. 
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Table 3. Parameters of the fit of equation ( |l8| ) against the results obtained from the completeness simulations for de Vaucouleurs profiles at redshifts z = 0.5 
and z = 0.75. Shown are the best-fit parameters of the formula given in equation (lla) po,mo, a, and /3. The normalisation po provides a rough estimate 
of the height of the plateau the detection probability reaches. The break magnitude mo can be used as an approximation of the 50% completeness limit, if a 
higher maximum detection rate is reached. 
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Table 4. Parameters of the fit of equation ( |l8| ) against the results obtained from the completeness simulations for de Vaucouleurs profiles at thi'ee distinct 
redshifts z G {1-0, 1.25, 1.50}. Shown are the best-fit parameters of the forinula given in equation ( p^ po, mo, a, and /3. The norinalisation po provides 
a rough estimate of the height of the plateau the detection probability reaches. The break magnitude mo can be used as an approximation of the 50% 
completeness limit, if a higher maximum detection rate is reached. The detection probability in field S6F1 at redshift z = 1.50 is too low to allow a reliable 
fit to the data. 



surveys as long as a threshold based detection algorithm and an 
adaptive isophotal magnitude measurement are used. 

We find a strong bias against bright profiles at higher 
redshifts in the rate of detection as well in the photometry of the 
objects. At low redshifts the completeness fraction is comparable 
to point-like sources. The measured magnitudes underestimate the 
true magnitudes at the bright end of the distribution. The sharp 
maximum in the centre of the profile results in an underestimate 
of the objects' radius leading to a too small isophotal radius in the 
magnitude measurement. At high redshifts the detection probability 
does not reach one even for the brightest objects, forming a plateau 
at lower values. This stems from the correlation of surface bright- 
ness and effective radius along the fundamental-plane relation. To- 
ward the bright end the object's surface brightness decreases, al- 
though the total magnitude increases. 

The detection of spiral galaxies at high redshift shows simi- 
lar deficiencies. Here the effective radius increases with magnitude, 
but the intrinsic mean surface brightness stays constant as predicted 
by the Freeman law. The rather flat light distribution of the expo- 
nential profile helps to estimate the object's radius correctly and 
accordingly to measure the magnitudes without large errors. 

These results are in gener al agreement with the pred iction of 
the visibility function theory ( Disne^^^^hilli^ ps 1983 ) and the 
brightness defects predicted by palcantoir7T9'98[ ) - except in the 
two most distant redshift-bins, where the predictions break down - 
using calculations of the lost-light fraction, as long as the effects of 
seeing are taken into account. 

With equations ( |lo| ) and (ju]) we provide a rather simple way 



to estimate the ~ 50% completeness limit of point-like sources 
represented by Moffat profiles. 
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S4F1 


0.95 


19.34 


0.05 


81.52 


0.89 


19.22 


0.10 


75.30 


S5F1 


0.99 


19.35 


-0.05 


101.58 


0.96 


19.26 


-0.16 


90.01 


S5F5 


0.99 


19.38 


-0.11 


107.76 


0.94 


19.31 


0.00 


95.69 


S6F1 


0.98 


18.98 


-0.21 


104.80 


0.90 


18.93 


-0.07 


86.86 


S6F5 


0.99 


19.56 


0.00 


113.58 


0.96 


19.48 


-0.03 


99.99 


S7F5 


0.98 


19.37 


-0.05 


80.93 


0.95 


19.27 


-0.11 


67.06 



Table 5. Parameters of the fit of equation ( |l8| ) against the results obtained from the completeness simulations for exponential profiles at at redshifts z = 0.5 
and 2 = 0.75. Shown are the best-fit parameters of the formula given in equation ( |l^ poi "^Oi ct. ™d /3- The normalisation po provides a rough estimate 
of the height of the plateau the detection probability reaches. The break magnitude mo can be used as an approximation of the 50% completeness limit, if a 
higher maximum detection rate is reached. 



Field 




z = 


1.0 






2 = 


1.25 






2 = 


1.5 






PO 


mo 


a 




PO 


mo 


a 




Po 


mo 


a 


/9 


S2F1 


0.84 


19.00 


0.32 


89.53 


0.69 


18.98 


1.38 


88.07 


0.55 


18.92 


3.27 


83.59 


S2F5 


0.86 


19.28 


0.41 


89.29 


0.80 


19.14 


0.24 


73.32 


0.67 


19.08 


0.48 


75.83 


S3F5 


0.92 


19.39 


-0.18 


90.77 


0.82 


19.31 


0.29 


81.43 


0.64 


19.27 


1.90 


83.07 


S4F1 


0.84 


19.14 


0.08 


61.14 


0.67 


19.08 


1.28 


65.21 


0.47 


19.12 


3.43 


72.29 


S5F1 


0.87 


19.20 


0.22 


95.03 


0.77 


19.14 


0.50 


79.92 


0.57 


19.11 


2.67 


82.73 


S5F5 


0.88 


19.22 


0.10 


87.80 


0.77 


19.17 


0.31 


81.99 


0.61 


19.13 


1.30 


84.03 


S6F1 


0.81 


18.84 


0.07 


82.68 


0.60 


18.86 


2.46 


93.53 


0.36 


18.90 


7.54 


120.92 


S6F5 


0.93 


19.42 


-0.07 


90.57 


0.85 


19.37 


-0.06 


88.12 


0.68 


19.32 


1.25 


93.77 


S7F5 


0.86 


19.26 


0.19 


73.53 


0.74 


19.22 


0.84 


75.93 


0.52 


19.24 


3.60 


72.31 



Table 6. Parameters of the fit of equation ( |l8[ ) against the results obtained from the completeness simulations for exponential profiles at three distinct redshifts 
2 £ {1-0, 1.25, 1.50}. Shown are the best-fit parameters of the formula given in equation (hsk po,mo, a, and /3. The normalisation po provides a rough 
estimate of the height of the plateau the detection probability reaches. The break magnitude mo can be used as an approximation of the 50% completeness 
limit, if a higher maximum detection rate is reached. 
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